Astronomy & Astrophysics manuscript no. al 


© ESO 2013 


March 4, 2013 





Redshift-space correlation functions in large galaxy cluster surveys 

P. Valageas 1 and N. Clerc 2 

1 Institut de Physique Theorique, CEA Saclay, 91191 Gif-sur-Yvette, France 

2 Max-Planck-Institut fur extraterrestrische Physik (MPE), Giessenbachstrasse 1, 85748 Garching, Germany 
Received / Accepted 



o 

(N 
c3 



(N 
(N 

o 

u 

6 



> 

oo 

in 
o 

(N 



X 



ABSTRACT 

Context. Large ongoing and upcoming galaxy cluster surveys in the optical, X-ray and millimetric wavelengths will 
provide rich samples of galaxy clusters at unprecedented depths. One key observable for constraining cosmological 
models is the correlation function of these objects, measured through their spectroscopic redshift. 

Aims. We study the redshift-space correlation functions of clusters of galaxies, averaged over finite redshift intervals, 
and their covariance matrices. Expanding as usual the angular anisotropy of the redshift-space correlation on Legendre 
polynomials, we consider the redshift-space distortions of the monopole as well as the next two multipoles, 21 = 2 and 
4. 

Methods. Taking into account the Kaiser effect, we develop an analytical formalism to obtain explicit expressions of 
all contributions to these mean correlations and covariance matrices. We include both shot-noise and sample-variance 
effects, as well as Gaussian and non-Gaussian contributions. 

Results. We obtain a reasonable agreement with numerical simulations for the mean correlations and covariance matrices 
on large scales (r > 10/i _1 Mpc). Redshift-space distortions amplify the monopole correlation by about 10 — 20%, 
depending on the halo mass, but the signal-to-noise ratio remains of the same order as for the real-space correlation. 
This distortion will be significant for surveys such as DES, Erosita and Euclid, which should also measure the quadrupole 
21 = 2. The third multipole, 21 — 4, may only be marginally detected by Euclid. 

Key words. Surveys - Galaxies: clusters: general - Cosmology: large-scale structure of Universe - Cosmology: observations 



1. Introduction 

The growth of structure in our Universe results from 
a competition between gravitational attraction and the 
mean Hubble expansion. Studying recent structures 
provides understanding of both these phenomena along 
with supplementary insights into the initial conditions 
describing primordial density fluctuations, assumed to 
be very close to Gaussian- distributed. Analyses of mas- 
sively large galaxy surveys dCple et al.l 120051: iMples et all 
120081 lEisenstein et al.l 120111: iBlake et al l l2011bl). galaxy 
cluster surveys (Bohrin ger et al.l 120011: iBurenin et al.l 



1972 



iPlanck Collaboration et all 120111 : iReichardt et al.l 
2012). Observed clusters are much more rare than individ- 



ual galaxies (current cluster surveys amount up to a few 
tens of objects per deg 2 ) but the link to their parent dark 
matter halo mass is better understood. Most of cluster cos- 
mological studies t hus rely on the abundance of clusters 



and its evolution (Borgani et al.l 120 01: Roz o et al 
iMantz et all 120101: ISehgal et al.l 120111: Icierc et al. " 



.20101: 
2012). 



However, very large o ngoing surveys ( e.g. Pierre et al.l201l[) 



20071: IKoester et all l2007al: IVanderlinde et al.l 120101: 



Mehrtens et al.l 1201 lh. weak-lensing studies (e.g. 
Benjamin et al.l 120071: iMassey et al.l [20071: iKilbinger etatl 
20091 : ISchrabback et all I2010T ) . and measurements of the 



or pl anned surveys (jPredehl et al . 2009; iRefregier et al.l 
I2010h promise to overcome the statistical limit and bet- 
ter constrain cosmological models by stuyding the spatial 
clustering of clusters over large volumes. 



baryonic acoustic oscillation peak (|Eisenstein et al.l 120051 : 
IBlake et"al"1 l2011at ISlosar et al.l 1201 ll ). are currently under 
way or planned in order to probe the physics of large scales 
in multiple complementary ways. 

Clusters of galaxies are defined as galaxy concentra- 
tions but they can also be regarded as a specific class 
of objects in their own. In the hierarchical framework 
of structure formation, they are the largest and last 
nonlinear objects being formed in a quasi-equilibrium 
state. Observationally, they are well-defined and system- 
atically identified as l arge X-ray emi t ters, thanks t o thei r 
hot gas phase (e.g. ISarazi n| 119881: iPierre et al.l 1201 llh 
and /or isolat e d (red ) galaxy overdensities (jAbell Il958t 
IKoester et all l2007bft . and/or via their imprint on the 
Cosmological Background Radiation (jSunvaev &: Zeldovichl 



This work follows our previous study ■: Valageas et al.l 

1201 lL hereafter Paper I) that analytically quantified covari- 
ance matrices for correlation functions in the framework of 
large galaxy cluster surveys. These calculations are essen- 
tial in order to account for uncertainties in clustering anal- 
yses. Most generally, these uncertainties combine a shot- 
noise term (arising from the discreteness of the objects) 
and a sample variance term (due to the finite volume of the 
Universe being surveyed), as well as mixed contributions. 
All are related to the survey design and we applied our for- 
malism to predict signal to noise ratios for observable quan- 
tities given different survey strategies (jPierre et alJ [2011: 
IValageas et al.l [20TTh . We presented our method in a way 
as close as possible to observational analyses, e.g. refering to 
the widely-used Landy-Szalay estimator for the two-point 
correlation function, integrated within fixed distance bins, 
averaging over all accessible redshift values, etc. Our re- 
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suits were found to be in good agreement with large N-body 
simulations. However, only real-space comoving coordinates 
were considered, therefore ignoring the existence of halo pe- 
culiar velocities. Because their radial component impacts 
object distance measurements via an additional spectral 
redshift, these peculiar velocities distort correlation func- 
tion estimates. 

The goal of the present paper is to estimate the correc- 
tions induced by such redshift-space distortions, still con- 
sidering galaxy clusters only, and to predict if such effects 
could be observed in present and future cluster surveys. 
Indeed, redshift-space distortions measurements present a 
growing interest as a cosmological probe. Schematically, 
this can be understood as follows: as peculiar velocities 
are linked to the gravitational well created by surround- 
ing matter, bulk flows in the large scale structure are di- 
rectly linked to the growth rate of structures. More pre- 
cisely, redshift-space distortions are sensitive to D + (a), the 
linear growth factor as a function of cosmic scale-factor a, 
and to / = dln_D + /dlna. In particular, observational re- 
sults have been obtained by va rious galaxy spectros copic 
surveys such as the 2dFGRS (iPeacock et al M200lh . the 



SDSS dCabre fc Gaztanagal [200g] Reid et al. 112012ft . 2dF- 
SDSS LRG and Q SO survey lRoss et al.1 l2007fi~~and the 
VVDS IGuzzo et all rf2008ft . 

In the following, we assume galaxy clusters to have been 
detected individually in a dedicated survey (e.g. in the X- 
ray band) and their redshift to have been subsequently de- 
termined, e.g. by measurements of their member galaxy 
spectra and averaging over the corresponding galaxy red- 
shifts. This is the most common - but observationally ex- 
pensive - method to determine cluster redshifts. In practi- 
cal applications, only a subsample of all cluster members 
can be observed due to observational constraints (limiting 
magnitude, compactness of the cluster in the instrumental 
focal plane, finite wavelength range of the spectrograph, 
identification of spectral lines, elimination of fore- and back- 
ground objects, etc.). This translates i nto statistical errors 
on the cluster redshif t estimate (e.g. iKatgert et al.l Il996t 
lAdami fc Ulmerl [20001. Most generally, spectroscopic red- 
shifts of galaxy clusters are considered to be precise to 
a(z) ~ 0.01(1 + z) or less. 

For the purpose of disentangling measurement errors 
from intrinsic covariances, we assume cluster redshifts to be 
known with infinite accuracy and neglect any kind of bias. 
Adding such sources of noise (that are not correlated with 
the shot-noise and sample-variance terms that we consider 
here) is not difficult as one only needs to add the relevant 
covariance matrices. We do not consider this second step 
here, which depends on the details of the instrument. 

As in Paper I, we consider in our calculations both the 
effects of shot-noise and sample variance. We expand the 
angular anisotropy of the redshift-space correlation func- 
tion on Legendre polynomials. Our formalism then leads to 
expressions for the monopole and the two multipoles 21 = 2 
and 21 = 4. We provide explicit expressions for the ex- 
pected quantities and their covariance matrices, taking into 
account the Kaiser effect and Gaussian plus non-Gaussian 
contributions. Keeping in mind the application to galaxy 
cluster surveys, all quantities are integrated within the 
< z < 1 redshift range. Halos are defined by their mass 
M and we describe the population by using a mass func- 
tion and a bias model from previous works. The three- and 
four-point correlation function are estimated by means of 



a hierarchical model (|Peebledll980tlBernsteinlll994ft . where 
those quantities only depend on combinations of two-point 
correlation functions. Comparisons with numerical simula- 
tions are made to assess the reliability of our analytical 
calculations. 

Therefore, our main improvements over previous works 
consist in the inclusion of Gaussian and non-Gaussian con- 
tributions in the computation of covariance matrices, as 
well as in the derivation of observables close to those ex- 
pected in practical studies focused on galaxy clusters. 

This paper is organized as follows. In Sect. [5] we first 
briefly describe the analytic models that we use to esti- 
mate the means and covariance matrices of redshift-space 
halo correlation functions, as well as the numerical sim- 
ulations that we use to check our results. Then, focusing 
on the Landy & Szalay estimator, we study the means of 
the various multipoles of the redshift-space correlation in 
Sect. [3J and their covariance matrices in Sect. EJ We com- 
pare our analytical results with numerical simulations for 
the two lowest multipoles (the u 2£ = 0" monopole and the 
u 2£ = 2" quadrupole). We apply our formalism to several 
real surveys in Sect.jSJ for the three multipoles generated by 
the Kaiser effect (21 = 0,2, and 4). We conclude in Sect. [5] 

We give details of our calculations in Appendices [X] and 
[Bl and we show the auto- and cross-correlation matrices of 
the "2£ = 2" multipole in Appendix [Cj 



2. Redshift-space halo density field 

Before we describe our analysis of the covariance matri- 
ces for halo number counts and correlation functions, we 
present in this section the analytic models that we use for 
the underlying halo distributions (mass function and bias, 
etc.) and the numerical simulations that we use to validate 
our results. 



2.1. Analytic models 

To be consistent with the numerical simulations, in Sects. [3] 
and [4] where we develop our formalism and compare our 
results with simulat ions, we use the WMAP3 cosmology 
(|Spereel et al.l l2007h . that is, Q m = 0.24, n dc = 0.76, 
n b = 0.042, h = 0.73, cr 8 = 0.77, n s = 0.958, and 
u>dc = — 1. In Sect. [5J where we apply our formalism to 
obtain forecasts for current and f uture surveys, we us e 
the more recent WMAP7 cosmology ijKomatsu et al.lfeoilft . 
that is, tt m = 0.274, fl dc = 0.726, fl b = 0.046, h = 0.702, 
erg = 0.816, n s — 0.968, and u>d c = — 1- 

In this paper, keeping in mind the study of X-ray clus- 
ters, we consider the number counts and correlations of 
dark matter halos defined by the nonlinear density con- 
trast S = 200. These halos are fully characterized by their 
mass, and we do not investigate the relationship between 
this mass and cluster properties such as the gas tempera- 
ture and X-ray luminosity. These scaling laws can be added 
to our formalism to derive the cluster number counts and 
correlations, depending on the quantities that are actually 
measured, but we keep a more general setting in this paper. 

2.1.1. Halo mass function and bias 



We us e the halo mass function, dn/d ln M, of Tinker et all 
), and the halo bias, b(M), of ITinker et al.1 (|2010h . 
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Then, we write the two-point real-space correlation function 
between two halos labeled V and "j" as 



6u = bibjtdxi -Xj|;z), 



(1) 



where £ is the real-space matter density correlation. This 
corresponds to the halo power spectrum 



P£jk;z) = b i b j P(k;z), 



(2) 



where P(k) is the matter density power spectrum. Here 
we used the approximation of scale-independent halo bias, 
which is valid to better than 10% o n scales 20 < x < 130/i _1 
Mpc (jManera fc Gazta naga 20l3), with a small feature on 
the ba ryon acoustic scale (r ~ 100/i _1 Mpc) of amplitude 
of 5% (jDesiacaues et al.ll201Q[ ). 

To simplify notations, we define the mean cumulative 
number density n of objects observed at a given redshift, 
within the mass bin [M_ , M + ] , as 



n(z) 



dn 



+ dM 

M dlnM 



(M,z), 



and the mean bias b as 



r M + A yf An 

b( : i 77( ; ) = / 



(3) 



(4) 



M- 



In the following we omit to write the explicit boundaries on 
mass since we only consider a single mass bin. It is straight- 
forward to extend all our results to the case of several mass 
bins, by writing the relevant mass boundaries and replacing 

factors such as b n 2 by bifLibjUj. 



2.1.2. Redshift-space correlation functions 

We denote by x and s real-space and redshift-space coordi- 
nates. They only differ by their radial component, because 
of the peculiar velocity along the line-of-sight, 



s = x 



aH 



(5) 



Here, e z is the unit vector along the line of sight and v the 
peculiar velocity. At linear order, the peculiar velocity field 
and the density contrast are related in Fourier space as 



v(k) = i a- 



dlnDi k ~ 



dt k 2 



<5(k), 



(6) 



where D + is the linear growing mode. From the conserva- 
tion of matter, (1 + 5' s ^)ds = (1 + <5) dx, one obta ins the 
redshift-space matter power spectrum (|Kaiserlll987t ) 



F( s )(k;z) = (l + //4) 2 P(fc;z), 



where 

d\nD A 



f 



d In a 



(7) 



(8) 



Here and in the following, we use a superscript (s) for 
redshift-space quantities. Following common practice, we 
assume that the halo velocity field is not biased. This yields 
at the same order the redshift-space halo power spectrum 



P h ^(k;z) = b\l + ^l) 2 P(k;z), 



(9) 



with j3 — f /b, where b is the mean bias of the ha lo popula- 
tion as defined in Eq.g]) (jPercival fc Whitdl2009T) . We write 
the associated redshift-space halo correlation function as 

Z h{s Hs;z) = b 2 ^(s;z), (10) 
where £( s ) (s) can be expanded as ()Hamiltonl 1 1 9921) 

2 



^)( S )=£^)( S )P 2 ,( M| 



e, • s 



(11) 



i=0 



where we introduced the functions 

^ 0) (s) =( 1 + f + T)f f A2(fc) Jo(fcs) ' (12) 



£ (s;4) (*) = 



4p 
3 



dk 



A 2 (k)j 2 (ks) 



8f3 2 f°° dk a2 „, . ., . 



(13) 



(14) 



where ji are the spherical Bessel functions and we intro- 
duced the logarithmic power A 2 (k) — 47rfc 3 P(fc). 

The function and its multipoles £( s;2£ ) depend on 
the halo population through the factor j3 (i.e., the bias can- 
not be fully factored out as in the real-space correlation 
([T])). However, it is convenient to introduce the auxiliary 
function £( s ) as in Eq. (fTU|) so that most expressions ob- 
tained in the following have the same form as in real space 
(which is recovered by setting j3 — 0, that is, £( s ) = £). 

Independently of the linear-order approximation 
only even multipoles can appear in the expansion (1111) be- 
cause of parity symmetry in redshift-space. Nonlinear cor- 
rections may add higher-order even multipoles that could 
be handled by a straightforward generalization of the ap- 
proach developed in this paper. 

The redshift-spac e distortions @ and (fTTj) only include 
the "Kaiser effect" (jKaiserl fl987l ). associated with large- 
scale flows (i.e., the collapse of overdense regions and the ex- 
pansion of underdense regions), which amplifies the power 
because the density and velocity fields are closely related 
as in Eq.Q. However, on smaller scales within virialized 
regions, there appears a large velocity dispersion that adds 
a decorrelated random component to t he relatio n ([5]). This 
yields the so-called "fingers of god" (jJacksonl Il972l ) that 
can be seen on galaxy surveys and give a damping factor to 
their redshift-space power spectrum. In this paper we focus 
on cluster surveys where this effect does not arise (as can 
be ch ecked in numerical simulations, ()Nishimichi fc Taruval 
l201l[ )l Indeed, because clusters are the largest nonlin- 
ear objects cluster pairs are still in the process of moving 
closer (or further appart in fewer cases) and have not had 
time to fall within a common potential well and develop a 
large velocity dispersion through several orbital revolutions. 
Therefore, we only need consider the Kaiser effect as in 
Eqs.Q and (fTTj) . This provides a good agreement with nu- 
merical simulations for s > 20/i _1 Mpc (within 10%, except 
for the multi pole 21 = 4 that show s larger deviations from 
linear theorv lReid fc White! ((201 lh . but we will see that this 
multipole is not very important for cluster studies), and for 
k < 0.1/iMpc" 1 in Fourier space (jPercival fc White! 120091 
iNishimichi fc Taruval[20Tl . 
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We assume that these expressions still provide a good 
approximation up to the weakly nonlinear scales that are 
relevant for studies of cluster correlations (Jj 6ft. _1 Mpc) and 
for numerical purposes w e use the popular fit ting formula to 
numerical simulations of ISmith et al.l (|2003l ) for the matter 
power spectrum. 

Nevertheless, the formalism that we develop in the 
following sections is fully general and does not rely on 
Eqs.(HU)-([Il]) and the truncation of the series (JTTJ) at 21 = 4. 
It can be applied to any model for the redshift-space corre- 
lation £( s )(s). If this leads to higher-order multipoles 6 S;2 ^ 
one simply needs to extend the sums to all required multi- 
poles in the various expressions given in Appendix [X] 

2.1.3. Three-point and four-point halo correlations 

The covariance matrices of estimators of the halo two- 
point correlation involve the three-point and four-point 
halo correlations. Therefore, we need a model f o r thes e 
higher-or der correla ti ons. A s in | Valageas et al.| (|201l[) , 
following iBernsteinl (119941) . ISzapudi fc Colombil Jl996), 
iMeiksin fc White! (11999^ " we use a hierarchical ansatz 
(jGroth fc Peebles! Il977t iPeeblesI Il980h and we write the 
real-space halo three-point correlation ( h as 

g 

Ci,2,3 = & i & 2&3 y [6,26,3 + 6,16,3 + 6,16,2] , (15) 

where we sum over all three possible configurations over the 
three halos labeled "1", "2", and "3", and the real-space 
halo four-point correlation if 1 as 



16 

6,26,36,4 



[6,26,36,4 
+ 11 eye], 



3 eye. 



(16) 



where "3 eye." and "11 eye." stand for three and eleven 
terms that are obtained from the previous one by permu- 
tations over the labe ls "1,2,3,4" of the four halos. As in 
IValaeeas etall (|201lh . we take for the normalization factors 
^3 and 54 their l arge-scale limit, which is obtained by per- 
turbat ion theory (|Peebleslll980tlFrvlll984t[Bernardeau et al.l 
[20023, 



S 3 = 



s 4 = 



34 



(n + 3), 



60712 62. . 7. . 2 
l323--Y ( " + 3) + 3 ( " + 3) ' 



(17) 



(18) 



where n is the slope of the linear power spectrum at the 
scale of interest. The model (fl~5)) - ([T5)) is the simplest model 
that is in qualitative agreement with large-scale theoretical 
predictions (which also give the scalings £ ~ £ 2 and rj ~ £ 3 , 
but where the normalization factors S3 and S4 show an 
additio nal angular dependen ce) and small-scale numerical 
results (jColombi et al.|[l996]) (which however show an addi- 
tional moderate scale-depende nce of the normalizatio n fac- 
tors S3 and Si). As shown in IValageas et al.l (|2011[ ). this 
model (|15p - (fT5|) provides reasonably good predictions for 
the covariance matrices of estimators of the cluster two- 
point correlations (which focus on weakly nonlinear scales). 

Next, for the redshift-space correlations £ h ( s ) and r] h ^ 
we use the same hierarchical ansatz (ll~5|) - (fT6l) . where in the 
right-hand-side we replace £ by the redshift-space two-point 
correlation £( s ) given in Eq. flTlT) . This redshift-space model 



has not been checked in details against numerical simula- 
tions in previous works and we compare our predictions 
with simulations in Sect. 2] below. We will check that it 
yields a reasonable agreement for covariance matrices (espe- 
cially along the diagonal) as for our purposes we only need 
a simple and efficient approximation. If one aims at pre- 
dicting the three-point and four-point redshift-space corre- 
lations for their own sake, up to a high accuracy, one should 
perform dedicated tests for these quantities and probably 
build more sophisticated models, but this is beyond the 
scope of this paper. 

This approximation allows us to exten d in a straight- 
forwar d fashion the formalism developed in Val ageas et al.l 
(|201lf ) to redshift-space estimators. The only difference 
comes from the angular dependence of the two-point cor- 
relation (TTTI) . that is, from the second and fourth spherical 

harmonics 6 s;2 H s )^2(m) an d 6 S;4 H S )-P4 (/■*)• I n particular, 
as clearly seen in Eq.©, by setting (3 — 0, which only leaves 
the monopole c ontribution £ = 0, we recover the real-space 
results given in IValageas et al.l f|201 lh . 

2.2. Numerical simulations 

To check our analytical model we compare our predic- 
tions with the num erical simulations that we used in 
IValageas et al.l (|201ll ) for real-space statistics. While we 
considered only real-space coordinates in our previous 
study, we make use of the peculiar velocity assigned to each 
halo in order to conduct our comparison in redshift-space 
coordinates. 

Here we briefly recall the main characteristics of 
these simulations. We use the high-resolut ion, all-sky, 
Horizon simulation run (jTevssier et al.l l2009h . It consists 
in a (2fa~ 1 Gpc) 3 N-body simulation carried out using the 
RAMSES code (|Tevssierll200"2h . The effect of baryons is ne- 
glected and the total number of particles in the simulation 
sets the total mass of the simulation. The mass resolution, 
which is the mass of each of the 4096 3 simulated particules, 
is ~ 10 10 MfTy Halos ar e found using the HOP algorithm 
(jEisenstein fc Hutlll998D . The (real-space) comovmg posi- 
tion and peculiar velocity of each halo are obtained by aver- 
aging over all particles making up its total mass. Comoving 
coordinates are further transformed into sky positions (two 
angles plus redshift) and redshift-space coordinates taking 
into account their radial velocity. As to ensure completeness 
of the simulation towards all directions, we only consider 
halos at z < 0.8. 

Simulated surveys are extracted from the full-sky light- 
cone catalogue by defining non-overlapping rectangular re- 
gions in angular coordinates. Each field covers an area of 
50 deg 2 . In order to minimize correlations between neigh- 
boring fields, we require a separation of 20 deg between two 
consecutive surveys, which results in a total of 34 such in- 
dependent fields distributed all across the sky. Only for the 
results discussed in Section 14.31 we also use a 10 deg gap 
(resulting in 138 fields). For the purpose of deriving the 
correlation function, auxiliary random fields are created by 
shuffling the angular coordinates of halos in the original 
data fields. This operation preserves the original redshift 
and mass distribution of the simulation. The density of ob- 
jects in the random fields is increased up to one hundred 
times that of original data fields to avoid introducing spu- 
rious noise. The Landy-Szalay estimator and its generaliza- 
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tion (see Section. 13.11 Eq. I2U1) provide an estimate of the 
two- point correlation and its multipoles in redshift space. 
As in I Valageas et al.l (|2011[) . the mean and covariance of all 
derived quantities are finally estimated by sample averaging 
over the extracted fields. 

3. Mean redshift-space two-point correlation 
functions 

Since we have in mind the application to cluster surveys 
and, more generally, to deep surveys of rare objects, we 
consider 3D correlation functions averaged over a wide red- 
shift bin (in order to accumulate a large enough number of 
objects), rather than the usual local 3D correlation func- 
tions at a given redshift. This means that the quantities 
that we consider, while being truly 3D correlations and not 
2D angular correlations, nevertheless involve integrations 
along the line of sight within a finite redshi ft interval. We 
follow the formalism and the notations of IValaeeas et all 
(|2011|) to derive the means and the covariance matrices of 
estimators of these integrated redshift-space correlations. 

3.1. Analytical results 

We focus o n the Landy fc Szalay es timator of the two-point 
correlation lLandv fc Szalavl (|l993r i. given by 



DD - 2DR + RR 
RR ' 



(19) 



where D represents the data field and R an independent 
Poisson distribution, both with the same mean density. In 
practice, before appropriate rescaling, the mean number 
density of the Poisson process R is actually taken as much 
higher than the observed one, so that the contribution from 
fluctuations of the denominator RR to the noise of £ can 
be ignored. The advantage of the form (|T9")) is that one au- 
tomatically includes the geometry of the survey (including 
boundary effects, cuts, etc.), because the auxiliary field R 
is drawn on the same geometry. In practice, our general- 
ization of the estimator (fT§|) to redshift-space multipoles 
writes as 



41, + 1 



^2 P 2^ (Mdd) 



(hi 



2}_^P 2U (n dr ) 

dr 



P2£ * 



(20) 



where each sum counts all the pairs, dd, dr, or rr, within 
the samples D and R, that are separated by a redshift-space 
distance s within the bin < s < Ri.+. As compared 
with the usual monopole estimator (|19l) , to each pair in the 
numerator we add the geometrical weight (4^ + l)P 2 £i (/■<) 
while in the denominator we keep the unit weight. Of 
course, for 4 = we recover Ea. ((T9"|) . 

Within our framework, we write this estimator t; R 2P , 
for the multipole 2£i of the mean correlation on scales de- 
limited by Ri t - and Ri,+ , integrated over some redshift 
range and mass interval, as 



dft 



dz J (An) 

AM 1 dh dn 



dM 
M~ 



ds'(44 + l)P2^') 



2 f. dx^ 2 



M> dlnMdlnM' 



dn 

(An) 



dM 
M~ 



x/ds / (4^ + l)P 2 ^( M / ) 



dM' dn 



dn 



M> dlnMdlnM 



77 + ^,0,(21) 



with 



Qi 



dz^V* 
dz 

dn 



dn 



(An) 



dM 
~M 



ds' 



dM' 



An 



dlnMdlnM' 



(22) 



Here we denoted j i ds' as the integral f R ' + ds' over the 3D 

spherical shell i, of inner and outer radii Ri t - and i2^+ in 
redshift-space, [i! = (e z ■ s')/s' is the cosine of the angle 
between the line of sight and the halo pair, x( z ) an d T3(z) 
are the comoving radial and angular distances, di^o is the 
Kronecker symbol (which is unity if £i — and zero oth- 
erwise), and dn/dlnM is the observed density of objects. 
Here and in the following, we note observed quantities by a 
hat (i.e., in one realization of the sky) to distinguish them 
from mean quantities, such as the mean comoving number 
density dn/dlnM. The two subscripts Ri and 2ii label the 
radial bin and the multipole. If we consider Nr radial bins 
and Ne multipoles there are Nr x Ni values of the index 
i (i.e., two different indices i and j may correspond to dif- 
ferent radial bins and multipoles, or to the same radial bin 
but with two different multipoles, and vice- versa). 

In Eqs. (f2"Tj) and (|2"2")l we considered a survey with a sin- 
gle window of angular area (An) and a single mass bin. 
It is straigtforward to generalize to several distant (uncor- 
rected) angular windows and to several mass bins. The 
redshift interval Az is not necessarily small, and to in- 
crease the statistics we can choose the whole redshift range 
of the survey. It is again straigtforward to generalize to 
several nonoverlapping redshift intervals, if they are large 
enough to neglect cross- correlations between different bins 
(Az ^ 0.1, see Fig. 11 in lValageas et al.l (|2011[ )h 

As in Eq. (|2"0")) . the counting method that underlies the 
first term in Eq. (|2"Tj) can be understood as follows. We 
span all objects in the "volume" (z, f2, InM), and count all 
neighbors at distance s', within the shell [-Ri,-, Ri t +], with 
a mass M' . We denote with unprimed letters the quanti- 
ties associated with the first object, (z, n, InM), and with 
primed letters the quantities associated with the neighbor 
of mass M' at distance s'. Thus, with obvious notations, 
dn/dlnM and dn/dlnM' are the observed number den- 
sities at the first and second (neighboring) points. In con- 
trast, the denominator Qi involves the mean number densi- 
ties dn/dlnM and dn/dlnM'. Therefore, Qi is not a ran- 
dom quantity so it shows no noise. Similarly, the difference 
between the terms associated with DD and DR is that in 
the former we have a product of two observed number den- 
sities, (dn/dlnM) x (dn/dlnM'), while in the latter we 
have a crossproduct between the observed and the mean 
number densities, (dn/dlnM) x (dn/dlnAf). 

Then, if £i ^ the second and third terms in Ea. (l2"Tl) 
vanish. This also means that for ii ^ the Landy & 
Szalay estimator (|20l) and the Peebles & Hauser estima- 
tor ((Peebles fc Hauserl Il974f) , which would read as £ = 
DD/RR - Si ii0 (with the implicit factor (4:£i + l)P 2 e z in 
DD), are equivalent (in the limit where the density of the 
auxiliary random sample R goes to infinity). 

The angular factor (4£i + l)P2i i (^') ensures that the 
estimator (|21l) extracts the multipole 2£i of the redshift- 
space correlation, as in (II 1 D . In our previous analysis of the 
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real-space correlation (Valageas ct al. 20lTJ) we only needed 
consider the monopole term, = 0, because higher multi- 
poles vanished. Here, because the redshift-space correlation 
is anisotropic higher multipoles are nonzero and contain 
valuable information. As seen in Eqs. (fT2"]) - (fT4")) , they con- 
strain the parameter j3 = f/bb, whence the growth rate of 
the density fluctuations and the bias of the halos. 
The quantity Qi also writes as3 



d X V 2 n 2 V l , 



where the volume V« of the z-shell is 

Air 

Vi(z) = y[i^ )+ (2) 3 - Ri,-(z)% 



(23) 



(24) 



which may depend on z. In practice, one would usually 
choose constant comoving shells, so that Vi does not de- 
pend on z. To obtain Eq. (j2"3"f we used that dn/dlnM and 
dn/dlnAf have no scale dependence (because they cor- 
respond to a uniform distribution of objects) and we ne- 
glected edge effects. (These finite- size effects are d i scuss ed 
and evaluated in Appendix B of Vala geas et al.l ()2011[ )h 
Next, the mean of the estimator reads as 



with 



d x v 2 b TrViGg^z), 



(4£ 4 + l)P 2£i ( M ')e (s) (s';^), 



(25) 



(26) 



which is the radial average of the 24-multipolc of ^ , over 
the 3D spherical shell associated with the radial bin i in 
redshift space. Indeed, from the multipole expansion 
this also reads as 



(s) _ dr 3r 2 , 2li) 



Hi.- _ 



(r). 



(27) 



3.2. Comparison with simulations 

We show in Fig. Q] the mean redshift-space and real-space 
halo correlations, (£r?q) an< ^ over ten distance bins. 

We find that because of the prefactor (1 + 2/3/3 + f3 2 /5) in 
Eq. (fT2"]) the "monopole" redshift-space correlation is greater 
than the real-space correlation by 10% for massive halos, 
M > W^h-Hlo, and by 18% for M > 2 x lO^/i^M©, in 
the redshift range < z < 0.8. The relative redshift-space 
distortion is smaller for more massive halos because they 
have a larger bias b. This decreases the factor /? = / fb. 

Although the results from numerical simulations are 
somewhat noisy they agree reasonably well with our analyt- 
ical predictions. They give the same order of magnitude as 
our predictions for both the redshift-space and real-space 
correlations. In particular, they fall within the expected 



1 We can neglect the impact of the change from x to s on 
background cosmological quantities such as the radial and an- 
gular distances \ and T>, which vary on scales on the order of 
c/Ho- Redshift-space distortions only have a significant effect on 
the correlation £' s ' of the halo population, which is studied on 
smaller scales ~ 20h~ 1 Mpc and becomes anisotropic. 



12 



10 



C 



0<z<0. 
l,=0 



M> 10 14 _h -'M r 



.M>2xl0 13 A-'M 




® 1 « 



10 



r [ h _1 Mpc ] 



10 s 



Fig. 1. The mean redshift-space (solid lines) and real- 
space (dashed lines) halo correlations, (Ir ) ano - (^i)- 
We consider ten comoving distance bins, within 5 < r < 
100/i _1 Mpc, equally spaced in log(r). We count all ha- 
los above the thresholds M* = 2 x 10 13 and 1O 14 /i _1 M , 
within the redshift interval < z < 0.8. We compare our 
analytical results with numerical simulations (crosses for 
the "monopole" redshift-space correlation and circles for 
the real-space correlation). The 3 — a error bars (slightly 
shifted to the left and to the right for the real-space and 
redshift-space correlations) are obtained from the analyti- 
cal covariance matrices derived in Sect. [4] 



3 — fx error-bars from our analytical predictions (these er- 
ror bars are obtained from the analytical covariance ma- 
trices derived in Sect. 0] below). This also shows that the 
accurate computation of such quantities remains a difficult 
task for numerical simulations (in this case it would re- 
quire many realizations of the simulated sky resulting in 
a heavy computational load) and it is useful to have an- 
alytical models for cross-checking. In particular, while the 

ratio (^r. o)/(^R>) ^ s greater than unity and smooth within 
the analytical model (it is actually scale-independent in 
our simple model (|12p) it shows spurious fluctuations in 
the simulation data for rare massive halos because of low 
statistics. 

For practical purposes, Fig. Q] shows that to constrain 
cosmology with the clustering of X-ray clusters it is neces- 
sary to take into account redshift-space distortions if one 
aims at an accuracy better than 10%. 

We show in Fig. [5] our results for the "2€ = 2" redshift- 

(s) 

space correlation £p 2 . This corresponds to the second 
term of the expansion ([lip over Legendre polynomials. As 
seen from Ea. (|13[) . this term is negative and its ampli- 
tude is larger for more massive halos because of the fac- 
tor t (4/3/3 + 4/3 2 /7) = 4/6/3 + 4/ 2 /7, which grows with 
the halo bias. However, the error bars being much larger 
for rare massive halos the associated signal is more difficult 
to measure. We obtain a good agreement with the simula- 
tions on large scales but on small scales below 10/i -1 Mpc 
our model seems to underestimate the amplitude of this 
quadrupole component of the redshift-space correlation. 
This may be partly due to higher-order contritbutions to 
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Fig. 2. The "2£ = 2 ,, redshift-space halo correlation (£^ 2 )- 
As in Fig.[TJ we consider ten comoving distance bins, within 
5 < r < 100/i _1 Mpc, equally spaced in log(r), and we 
count all halos above the thresholds M* — 2 x 10 13 and 
10 14 /i _1 M Q , within the redshift interval < z < 0.8. We 
compare our analytical results with numerical simulations 
(crosses). The 3 — a error bars (slightly shifted to the left 
and to the right for the real-space and redshift-space cor- 
relations) are obtained from the analytical covariance ma- 
trices derived in Sect. S) 



Eq. that modify the relationship between real-space and 
redshift-space power spectra. For instance, within pertur- 
bation theory next-to-leading terms of order P| arise and 
have been found to increase the amplitude of the redshift- 
space power and to improve the agreement with simulations 
(|Taruva et alJl2010t iNishimichi fc Taruvall201lh . However, 
we do not investigate such modifications in this paper. 



(s) 

We do not plot the last "2£ = 4" term, £V 4 , because our 
analytical computations show that it is very small (smaller 
than 0.1 on the scales shown in Figs. Q] and [5]) and within the 
error bars. However, we have checked that our model agrees 
with the simulations for r > lOh" 1 Mpc (on smaller scales 
our model again seems to underestimate this correlation). 



4. Covariance matrices 

4,1. Explicit expressions 



(s) 

We now consider the covariance matrices C) ' of the esti- 



mators ^jf , 



. , defined by 



a 



i{s) 



i{s) 



(28) 



To simplify notations, we do not write the labels Ri and 2£i 
in the covariance matrices and the subscript i again refers 
to both the radial bin, [R^Rjj , , and the multipole, 2£i. 

The results of I Valageas et all ( 201lh directly extend to 
redshift space and we obtain the expression 



with (see also lLandv fc Szalavl (|l993f l: ISzapudil (|200lh : 
iBernardeau et all (|2002j ) 



a 



«(2) 



2(4l i + l)(4l j + l) 
dMii dn dn 



dSli dMi 

(a^)~mT 



dsi. 



M/ dlnAf, dlnAf,-/ 



Ms) 



XP2^ (Hi>)P2ti{Hi>) 



c 



(s)(3) _ 4(44 + 1) (41, + 1) 
ds, d ^' 



dn 



2 dn, dMi 

dn dn 



ds,y - 



d In Mi d In M v d In My 



a 



(s)( 4) = (4£ t + l)(4£ J + l) 
dn dn 



dOi dMi 

(An) Mi 



ds,,, 



(30) 



M^ 



(31) 



dM e 
Mv 



d In M d In Mv 



, dn, dM, /■ dM ; - 



dn 



d In Afj d In Mj 



„Ms)Ms) Ms) 
si;j' sj'jj T" 'n,i';j,y 



(32) 



Here we use the same notation as in I Valageas et al.l (|201l| ). 
that we recalled below Eq. (l2"2"]) . Thus, labels i and j refer 
to objects that are at the center of the V, and Vj shells, 
whereas labels i' and j' refer to objects that are within the 
Vi and Vj shells. The pairs {i,i'} and {j,j'} are located 
at unrelated positions in the observational cone, which we 
integrate over. 



the Kronecker 



In Eg. (1501) we note 8n u R j 
with respect to the radial bins. Thus, 8R it R 
\Ri_-,Ri 



is 



symbol 
unity if 



[Rj-,Rj,+] and zero otherwise (we consider 
non-overlapping radial bins), independently of the multi- 
poles ii and tj (which may be equal or different). This 
shot-noise contribution arises from cases where the pairs 
{i,i'} and {j,j'} are the same (i.e., they involve the same 
two halos), which implies the same pair length. 

The label refers to quantities that involve n distinct 
objects. Thus, the contributions and arise from 
shot-noise effects (as is apparent through the prefactors 
l/(An)), associated with the discreteness of the number 
density distribution, and they would vanish for continuous 
distributions. However, they also involve the two-point and 
three-point correlations, and as such they couple discrete- 
ness effects with the underlying large-scale correlations of 
the population. In case of zero large-scale correlations, 
remains nonzero because of the unit factor in the brack- 
ets and becomes a purely shot-noise contribution, arising 
solely from discreteness effects. The contribution is a 
pure sample- variance contribution and does not depend on 
the discreteness of the number density distribution (hence 
there is no l/(An) prefactor). 

We can reorganize these various terms into three con- 
tributions, 



a 



c 



(a)(2) 
i,3 



c: 



W(3) 



00(4) 



(29) C$=C, 



-.(s)(G) 



c 



W(0 

i,3 



c 



i,3 ' 



(33) 
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where Cj*- G ^ gathers all terms that only depend on the 
two-point correlation, arises from the three-point 



correlation (i.e., the last term in Eo . (|3"Tj) ) . and C, 
arises from the four-point correlation (i.e., the last term 
in Eq.(j3U)). We obtain for the Gaussian contribution, 



(Afi)Q. 
4 



(Afl)Q i Q j 
2 



QiQj 



(34) 



which involves terms that are constant, linear and quadratic 
over fW. Here we introduced the mean 



with 



(35) 



dV 



(44+1)(4^ + 1)P 2 ^( M ')-P2^(m') 



x^)(s';z). (36) 
From the multipole expansion this also reads as 

/■«*,+ dr 3r 2 



41^,2*, = (4^ + l)(4^ + l) 



xE( 2 o f ^'T^W' 

£=0 ^ ' 



(37) 



which involves the square of the Wigner 3-j symbol. 

Using the hierarchical ansatz described in Sect. 12.1.31 
the contribution * associated with the three-point cor- 
relation writes as 



*' J (Afi)Q 4 Q_ 



|d X P 2 6 3 n 3 HV^ 



(38) 



while the contribution C) ■ associated with the four- 
point correlation writes as 



C 



0)W 



_i_oc( s )t A s ) , t( s ) c( s ) f( s ) < t( s ) A s ) A s ) (qn\ 

Here the quantities that are overlined are various geomet- 
rical averages of the redshift-space two-point correlation 
functions. For instance, the geometrical average introduced 
in the first integral in Eg . (|3"4")) writes as 



= (44+1x4^+1) l ^fPtzM f^i; 



x$ (s) (s r 



(40) 



We give the explicit expressions of these various terms in 
Appendix |21 and we describe how angular integrations can 



10 s 
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C 0.1 



io- 
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Fig. 3. The redshift-space (solid line, for 21 = 0) and 
real-space (dashed line) covariance matrices, C^'J and C^j, 
along the diagonal, i = j. We show the results obtained for 
halos in the redshift range < z < 0.8 with an angular 
window of 50 deg 2 . 



be performed using the spherical harmonic decomposition 
(fTTj) . as well as the flat-sky and Limber's approximations. 
All term s involve the redshift-space correlation except 
for £ cy i that only depends on the power along the trans- 
verse directions (and is not affected by redshift-space dis- 
tortions). 



4.2. Covariance matrix of the "monopole" redshift-space 
correlation 

(s) 

In this section we study the covariance matrix C\ ! of the 
"monopole" redshift-space correlation, £zj , that is, i{ = 



0. 



4.2.1. Comparison between redshift-space and real-space 
covariance matrices 

We show in Figs. [3] and H] the redshift-space and real-space 

(s) 

covariance matrices G-j and Cij, for halos in the redshift 
range < z < 0.8. The redshift-space covariance matrix is 
greater than the real-space covariance matrix by about 11% 
along the diagonal and by 20% for off-diagonal entries, for 
massive halos, M > lO^h^M®. For M > 2 x HF/i^Mq, 
these amplification factors are larger, 33% along the diago- 
nal and 40% for off-diagonal entries. This could be expected 
from the larger value of the redshift-space correlation func- 
tion noticed in Fig.Q] This is due to the amplification of the 
redshift-space power spectrum through the Kaiser effect, as 
shown by Eq.(j9]). Again, the relative effect is smaller for 
rare massive halos because of their larger bias b, which de- 
creases the associated factor (3 — f jb. The amplification of 
the covariance matrix is smaller along the diagonal, where 
the covariance matrix includes a pure shot-noise term (first 
term in Eg . AMI) ) that is independent of ^ s - ) and is not am- 
plified by the Kaiser effect. 

We obtain a reasonably good agreement with the nu- 
merical simulations for both redshift-space and real-space 
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Fig. 4. The redshift-space (solid line, for 21 — 0) and real- 
space (dashed line) covariance matrices, C^j and Cij, as 
a function of j, for i = 4 associated with the distance bin 
12.3 < r < 16.6/i _1 Mpc. We show the results obtained for 
halos in the redshift range < z < 0.8 with an angular 
window of 50 deg 2 . 



Fig. 6. The relative importance of the non-Gaussian con- 
tributions Cg (0 (solid line) and CjJ (n) (dashed line) to 

( s) 

the redshift-space covariance matrix ■ , along the diag- 
onal, i = j. We show the results obtained for halos in the 
redshift range < z < 0.8 with an angular window of 50 
deg 2 . 



10 



x > 
O 
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Fig. 5. The signal-to-noise ratio of the redshift-space (solid 
line) and real-space (dashed line) two-point correlations. 
We show the results obtained for halos in the redshift range 
< z < 0.8 with an angular window of 50 deg 2 . 



We show in Fig. [5] an estimate of the signal-to-noise ratio 
of the redshift-space and real-space two-point correlations, 
defined as ) / yJCTJ- A more precise analysis would re- 
quire using the full covariance matrix, but this ratio should 
allow us to check whether redshift distortions, which am- 
plify both the mean correlation and its covariance matrix, 
degrade or not the signal-to- noise ratio. Thus, Fig. [5] shows 
that the signal-to-noise ratio is not strongly affected by 
the Kaiser effect. In fact, the amplification of Cr o over- 



( s) 

weights the amplification of C\ J and the signal-to- noise ra- 
tio slightly increases. This can be understood from the fact 
that the contributions to c\ s ^ that are cubic over ^ (or of 
higher order) , which in our framework only come from the 
four-point correlation contribution (|39[) to the pure sample- 
variance term, are less important than the mixed shot-noise 
and sample- variance contributions that are constant or lin- 
ear over £( s ). 

The results found in Fig. [5] mean that redshift distor- 
tions do not prevent the measure of the clustering of X-ray 
clusters and its use to constrain cosmology, since the signal- 
to-noise ratio remains as large and even slightly greater. 



covariance matrices, for diagonal and off-diagonal entries. 
The simulation data are somewhat noisier for off diagonal 
entries and show spurious oscillations in Fig.|4]that are not 
physical but due to low statistics, especially for rare mas- 
sive halos. This shows that analytic models are competitive 
tools to estimate such covariance matrices. In particular, 
for data analysis purposes it can be helpful to have smooth 
covariance matrices, which make matrix inversion more re- 
liable. 

In any case, redshift-space distortions only have a mod- 
erate impact on the covariance matrices and do not change 
their magnitude. This means that most of the real-space 
results obtained in lValaeeas et al.l (|201 lh remain valid. 



4.2.2. Importance of non-Gaussian contributions 

We show in Figs. E and [7] the ratios C^) iC) /C^) and 

give the relative contribution to the 
covariance matrix of the terms associated with the three- 
point correlation function, Eq. (|38[) . or with the four-point 
correlation function, Ea. ((39|) . For low- mass halos the 4-pt 

(s)(ri) 

contribution C\ \J is larger than the 3-pt contribution 

Cj S j^\ whereas the ordering is reversed for high-mass ha- 
los. This agree s with the results obtained in real-space in 
IValaeeas et al] (|201lD (see their Figs. 15 and 16). This be- 
havior is due to the greater importance of shot-noise ef- 
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Fig. 7. The relative importance of the non-Gaussian con- Fig. 9. The covariance matrix C\ along the diagonal 



tributions C, 



W(0 



(solid line) and C^ (7?) 



(dashed line) to 



the redshift-space covariance matrix C\ , along one row 



(i = 4), associated with the distance bin 12.3 < r < 
16.6/i~ 1 Mpc. We show the results obtained for halos in the 
redshift range < z < 0.8 with an angular window of 50 
deg 2 . 



i — j, for the "21 = 2" mutipole of the redshift-space cor- 
relation. We show the results obtained for halos in the red- 
shift range < z < 0.8 with an angular window of 50 deg 2 . 
Points from numerical simulations are obtained from 138 
fields separated by a 10 degree gap (squares) or from 34 
fields separated by 20 degree gap (crosses). 



10 1 



fects for rare massive halos. This gives more weight to the 

3- pt contribution, which is part of the coupled shot-noise- 
sample- variance contribution C (s)(3) in Ea.pT]). than to the 

4- pt contribution, which is part of the pure sample- variance 
contribution C^ 4 > in Eq.©. For the same reason, these 
non- Gaussian contributions are relatively smaller along the 
diagonal for rare halos, as seen on the fourth bin j = i in 
Fig. [71 

Figs. |5]and [7] show that the non-Gaussian contributions 
can make up to 20 — 60% of the full covariance matrix 
as soon as one bin i or j corresponds to scales on the 
order of lOh^ 1 Mpc. Their relative contribution is also 
somewhat larger for off-diagonal terms. Therefore, although 
these terms do not change the order of magnitude of the 
covariance matrix, it is necessary to include them to obtain 
accurate estimates. 



4.2.3. Correlation matrices 

We show in Fig. \E\ the redshift-space correlation matrices 
defined as 



K {s) = 



C 



i,3 



(41) 



We plot the correlation matrices associated with the 
Gaussian part (f34|) and with the full matrix (j33)) . In agree- 
ment with Fig. [71 we can see that the non-Gaussian con- 
tributions, associated with the 3-pt and 4-pt correlations, 
make the correlation matrix less diagonal, especially for 
low-mass halos. As explained above, this is because rare 
massive halos show larger pure shot-noise effects that give 
an additional weight to the diagonal entries. In any case, 
we can see that it is important to take into account 
non-Gaussian contributions to obtain accurate estimates. 
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Fig. 10. The covariance matrix C-j as a function of j, 
for i = 4 associated with the distance bin 12.3 < r < 
16.6/i _1 Mpc, for the a 2i = 2" mutipole of the redshift- 
space correlation. We show the results obtained for halos 
in the redshift range < z < 0.8 with an angular window 
of 50 deg 2 . 



Moreover, off-diagonal terms cannot be neglected (espe- 
cially for low-mass halos) hence it is necessary to use the 

(s) 

full covariance matrix C) ' to draw constraints on cosmol- 
ogy from the analysis of X-ray clusters surveys. 



4.3. Covariance matrix of the '11 = 2" redshift-space 
correlation 

(s) 

We now study the covariance matrix C\ ■ of the 21 = 2" 



(s) 

redshift-space correlation, £ R 2 , that is, 21 



2L 



2. 
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Fig. 8. Contour plots for the redshift-space correlation matrices Tcfj , defined as in Eq . (|4"Tj) . There are ten distance 
bins, over 5 < r < 100/i _1 Mpc, equally spaced in log(r), as in previous figures. We consider halos in the redshift range 
< z < 0.8, within an angular window of 50 deg 2 , above the mass thresholds M > 2 x 10 13 /i _1 A/ Q in the upper row, 
and M > 10 14 /i~ 1 Af Q in the lower row. Left panels: correlation matrix TZ\ S ^ G ^ associated with the Gaussian part (|34p of 
the covariance matrix. Right panels: full correlation matrix, associated with the full matrix (|33[) . 



We compare in Figs. 151 andlTUlour analytical model with 
numerical simulations, for the covariance matrix along the 
diagonal and along one row. We consider two analysis of the 
numerical simulations, where we use either 138 nonoverlap- 
ping fields separated by a 10 degree gap (squares) or 34 
fields separated by a 20 degree gap (crosses). This gives 
a measure of the statistical error of the simulation data, 
apart from the systematic finite-resolution effects. Along 
the diagonal (Fig. [9]) our predictions show a reasonable 
agreement with the simulations, which appear reliable with 
low noise. Along one row (Fig. ITUf the simulations give 
higher off-diagonal entries than our model but they also 
show a significant statistical noise. In view of these uncer- 
tainties we obtain a reasonable agreement on scales above 
10/i _1 Mpc. On smaller scales, our model may underesti- 
mate the off-diagonal covariance because of nonlinear ef- 
fects, as suggested by the comparison for the mean cor- 
relation itself shown in Fig. [5J Therefore, on scales above 
10/i _1 Mpc analytical models such as the one presented in 
this paper can provide a competitive alternative to numer- 
ical simulations to estimate quantities such as covariance 
matrices that are difficult to obtain from numerical simula- 
tions (as they require large box sizes with a good resolution 
and time-consuming data analysis procedures). 



We show in Figs. [TT] and the relative contribution 
to the covariance matrix of the terms associated with the 
3-pt and 4-pt correlation functions. The comparison with 
Figs. |5]and [7] shows that these non-Gaussian contributions 
are relatively less important along the diagonal but more 
important on off-diagonal entries than for the monopolc 
u 2£ = 0" case. Thus, they make up to 40 - 60% of the full 
covariance matrix over a large range of scales. 

For completeness we show in Appendix[C]the correlation 
matrices of this "2£ = 2" redshift-space correlation, as well 
as the cross-correlation between the "2£ = 0" and u 2£ = 2" 
components. We find that the u 2£ — 2" covariance matrix 
is significantly more diagonal than the u 2£ = 0" covariance 
matrix and that the cross-correlation matrix between both 
components is quite small (~ 0.05). This means that the 
full covariance matrix in the {Ri, £{} space is approximately 
block-diagonal (we may neglect entries with £i ^ ij). Then, 

we can decouple the analysis of an d 2 - 

5. Applications to real survey cases 

As in IValaeeas et al.l (|201lD . we now study the correlation 
functions that can be measured in several cluster surveys. 
We consider three surveys on limited angular windows: 
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Fig. 11. The relative importance of the non-Gaussian con- 
tributions C^ )(C) (solid line) and C^ (7?) (dashed line) to 

(s) 

the redshift-space covariance matrix C\ ■ , along the diag- 
onal, i = j. We show the results obtained for halos in the 
redshift range < z < 0.8 with an angular window of 50 



.MO)) 



deg . (The contribution C %J 
not appear because it is smaller than 10~ 2 .) 



for M > 10 u /i- 1 M m does 



o 

"5 o.i 



10" 3 
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Fig. 12. The relative importance of the non-Gaussian con- 
tributions Cj)^ (solid line) and C^j' (dashed line) to 



(s) 

the redshift-space covariance matrix C\ ■ , along one row 
(i = 4), associated with the distance bin 12.3 < r < 
16.6/i _1 Mpc. We show the results obtained for halos in the 
redshift range < z < 0.8 with an angular window of 50 
deg 2 . 



- The XXL survey (jPierre et al.ll201lh is an XMM Very 
Large Programme specifically designed to constrain the 
equation of state of the dark energy by using clusters of 
galaxies. It consists of two 25 deg 2 areas and probes mas- 
sive clusters out to a redshift of ~ 2 We consider the "CI 
selection function" give n in FigJU of lValageas et all (pOll 
sec also lPacaud et all (|2006L 120071) 1 



Fig. 13. The mean "monopole" redshift-space (solid lines) 

and real-space (dashed lines) correlations, (i^ ) and 

(£r»)> over ten comoving distance bins within 5 < r < 
100/i _1 Mpc, equally spaced in log(r). We integrate over 
halos within the redshift interval < z < 1, for the XXL, 
DES, and SPT surveys. For DES we consider the mass 
thresholds M > 5 x 10 13 /i _1 M Q and M > 5 x W^H^Mq 
(smaller error bars), and for SPT the mass threshold M > 
5 x 1Q h Mq (larger error bars shifted to the right). The 
error bars show the diagonal part of the redshift-space co- 



variance f^(a) 



C 



1 5 



10 



a 

CO 



0<z< 1 



M>5xl0 14 ii-»M E) 

(Erosita & Euclid), 




1 O 



r [ h _1 Mpc ] 
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Fig. 14. The mean "monopole" correlations, ) (solid 

lines) and } (dashed lines) , as in Fig. [13] but for all- 
sky surveys. From top to bottom, we show a) halos above 
5 x lO 14 /i _1 M0 in either Erosita or Euclid, b) halos detected 
by Planck, c) halos above 5 x 10 13 h~ 1 M Q in Euclid, and d) 
halos detected by Erosita. 



- The Dark Energy Survey (DES) is an optical imaging 
survey to cover 5, 000 deg 2 with the Blanco 4- meter tele- 
scope at the Cerro Tololo Inter- American Observatorj0. We 
consider the expected mass threshold M > 5 x 10 13 H~ 1 Mq, 



https:/ /www. darkenergysurvey.org/index.shtml 
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shown in 



The ratio ((£«,) 

Fig.m 



10 s 



for the surveys 



are expected ove r the whole sky. W e assume the selec- 
tion function b y Melin et al.l ()2006f) . shown in Fig.J.l of 
IValaeeas etall (|2011[ ). 

- For Erosita, a simple flux limit is currently assumed 
as an average over the whole sky: 4 10~ 14 ergs -1 cm~ 2 
in the [0.5-2] keV band (|Predehl et al.ll2009ft. The associ- 
ated s election function is shown in Fig.J.l of Valageas et al.l 
1)201 lh . This would yield 71, 907 clusters out to z = 1. 

- For Euclid, we follow the prescription of the Euclid 
Science Book for the cluster optical selection function 
and adopt a fixed mass threshold of 5 x 1O 13 /i _1 M0 
(|Refregier et al.ll2010D . 

In the following the error bars shown in the figures are 
the statistical error bars associated with the covariance ma- 
trices studied in Sect. |4j (i.e., we do not include other sources 
of uncertainties such as observational noise) . In the case of 
full-sky surveys this "sample-variance" noise does not van- 
ish because even full-sky surveys up to a given redshift only 
cover a finite volume. In this limit the "sample variance" is 
usually called the "cosmic variance" . 
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Fig. 16. The ratio ((^o) - (£r< )) /<?£<.>) , for th e surveys 
shown in Fig. Q31 



io a 



as well as the subset of massive clusters M > 5 x 
10 u h- l M®. 

- The South Pole Telescope (SPT) operates at millime- 
ter wavelength^. It will cover some 2, 500 deg 2 at three 
frequencies, aiming at detecting clusters of galaxies from 
the Sunyaev-Zel'dovich (S-Z) effect. We consider a mass 
threshold of 5 x 10 14 /i _1 M s (jVanderlinde et alJl2010h . 

We also consider three all-sky surveys. In practice, the 
total angular area of such surveys is not really 4ir sterad 
since we must remove the galactic plane. In the following, 
for Planck we consider the two-sided cone of angle S = 75 
deg (i.e., \b\ > 15 deg), which yields a total area Ail ~ 
30576 deg 2 . For Erosita and Euclid we take 9 S = 59 deg 
(i.e., |6| > 31 deg), which corresponds to a total area that 
is about one-half of the full sky, AO, ~ 20000 deg 2 . 

- Planck operates at nine frequencies, enabling an effi- 
cient detection of the cluster S-Z signature but has a rather 
large PSF (5'-10'). Some 1625 massive clusters out to z = 1 



5.1. "Monopole" redshift-space correlation 

We show the redshift-space and real-space correlations 
(£r o) an d (ii) that can be measured with these surveys 
in the redshift range < z < 1 in Figs. [T3]and[T31 In each 
case, the redshift-space correlation is slightly higher than 



the real-space correlation. The error bars cr^s) 



C 



(<0 



http:/ /pole. uchicago.edu/ 



correspond to the diagonal entries of the redshift-space co- 
variance matrix. 

Redshift-space distortions can be seen on these figures 
but in several cases they are smaller than the error bars. To 
see more clearly the impact of these redshift-space distor- 
tions, we show the ratio ((£^ ) — (^Ri))/^*) in Figs. [TBI 

and [16] This gives the difference between the "monopole" 
redshift-space and real-space correlations in units of the 
standard deviation cr^s) , that is, of the error bar of the 
measure. 

Among limited-area surveys, the amplification of the 
cluster correlation function by the Kaiser effect will be 
clearly seen for halos above M > 5 x 10 13 ft- _1 Af Q in 
DES. The effect is negligible for more massive halos, M > 
5 x 1O 14 /i- 1 M , in DES and SPT, because of their larger 
error bars (see Fig. [T3l . For the "CI clusters" measured 
in the XXL survey redshift distortions are marginally rel- 
evant: they are on the order of 0.4ct^( S ) in distance bins 

around 15ft.~ 1 Mpc, so that by collecting the signal from all 
distance bins they would give a signal-to-noise ratio of order 
unity. 

Among full-sky surveys, the redshift-space amplification 
will have a clear impact for the full cluster populations 
measured in Erosita and Euclid, and a marginal impact 
for Planck. Again, subsamples of massive clusters, M > 5 x 
10 14 /i _1 M Q , lead to large error bars that hide this redshift- 
space effect. 

In practice one does not measure both real-space and 
redshift-space correlations from a cluster survey, because 
one does not have a map of the cluster velocities. Therefore, 
the redshift-space amplification of the spherically aver- 
aged correlation ) is degenerate with the halo bias 
b (which is difficult to predict with a high accuracy). To 
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Fig. 17. The mean "21 — 2" redshift-space correlation 
(£^2)- As in Fig. [T3] we consider ten comoving distance 
bins within 5 < r < 100/i _1 Mpc, equally spaced in log(r), 
and we integrate over halos within the redshift interval 
< z < 1. We show our results for halos above M > 
5 x 10 13 /i _1 Af© for DES. The error bars show the diagonal 



part of the redshift-space covariance a^a) = 



C 



-0.5 



1.5 



Erosita 




Euclid 

M>5x 10 13 h -'M„ 



1 O 



r [ h _1 Mpc ] 



10 z 



Fig. 18. The mean "2£ — 2" redshift-space correlation 
(|jj? 2 )) as ln Fig. [T7] but for all-sky surveys. We show our 
results for all halos detected by Erosita (upper curve) and 
for halos above 5 x 10 13 /i _1 Mq in Euclid (lower curve). 



measure redshift-space distortions one needs to measure 
the anisotropics of the redshift-space correlation, that is, 
the higher-order multipoles £ (s;2) and £ (s;4) of Eq.([TT]). This 
provides a distinctive signature that breaks the degeneracy 
with the halo bias and yields a further constraint on cos- 
mology through the factor j3, that is, the growth rate / 
of the density fluctuations defined in Eq.©. We consider 
these higher-order multipoles in the following sections. 



Fig. 19. The mean u 2£ — 4" redshift-space correlation 
4}- As in Fig. [T3l we consider ten comoving distance 
bins within 5 < r < 100/i _1 Mpc, equally spaced in log(r), 
and we integrate over halos within the redshift interval 
< z < 1. We show our results for halos above M > 
5 x lO 13 ^ 1 !/© for DES (left error bars) and Euclid (right 
error bars). 

5.2. "21 = 2" redshift-space correlation 

We show in Figs. [17] and QH the mean "2£ = 2" redshift- 
space correlation 2 }- We consider the same surveys as 
in Sect. I5.1l but we only plot our results for the cases where 
the signal-to-noise ratio is larger than unity. In agreement 
with the analysis of Sect. [31 this "quadrupole" term is more 
difficult to measure than the "monopole" term studied in 
Sect. 15. ll and only DES, Erosita, and Euclid can obtain a 
clear detection. Even for these surveys, the correlation can 
only be measured for the full halo population and error bars 
become too large if one restricts to the sub-sample of rare 
massive halos above 5 x lO 14 /i~ 1 M0. Nevertheless, it will be 
useful to measure this u 2£ = 2" redshift-space correlation in 
these surveys because this should tighten the constraints on 
cosmology. Indeed, a simultaneous analysis of (Cr^o) anu 

(oj. 2) constrains the factor j3 — f /b, whence the halo bias 
and the growth rate of density fluctuations. 



5.3. "21 = 4" redshift-space correlation 



We show in Fig. [19] the mean u 2£ = 4" redshift-space cor- 
The signal is even weaker than for the 
2" correlation and only DES and Euclid may be able 



relation (C^ 1 4 ) 
"2£ ■ 



to obtain a detection. DES is unlikely to provide an accu- 
rate measure but it should still give useful upper bounds. 
Euclid may obtain a significant measure of the overall am- 
plitude by gathering the information from all distance bins. 
This would provide an additional constraint on /3 to the one 
provided by the mutipole u 2£ = 2" of Fig. QB 



6. Conclusion 

In this paper we have generalized our previous study of 
the real-space correlation functions of clusters of galaxies 
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(|Valageas et al.ll20lTl ) to include redshift-space distortions. 
On large scales they lead to an anisotropic correlation func- 
tion because of the Kaiser effect, due to the correlation be- 
tween density and velocity fields. There are no "fingers of 
god" because clusters are the largest bound objects, as op- 
posed to galaxies that can have performed several orbital 
revolutions within a larger halo and show a large small- 
scale velocity dispersion. Then, taking into account this 
Kaiser effect at leading order and using a simple hierarchi- 
cal model for the three-point and four-point halo correla- 
tions, we have developped an analytical formalism to obtain 
explicit expressions for the mean redshift-space correlations 
and their covariance matrices. We include both shot-noise 
and sample-variance effects, as well as Gaussian and non- 
Gaussian contribution s. This is a direct ex tension of the 
formalism presented in Valagea s et al.l (|201ll) for real-space 
correlations. 

Expanding as usual the angular anisotropy of the 
redshift-space correlation on Legendre polynomials, we 
have considered the first three nonzero multipoles: the 
monopole 2£ = (i.e. the spherically-averaged correlation), 
the quadrupole, 2£ = 2, and the next multipole, 2£ = 4. 
Higher-order multipoles should arise from nonlinear cor- 
rections to the Kaiser effect. Although they are not stud- 
ied here, which would require a more complex nonlinear 
modeling of the Kaiser effect itself, their analysis could be 
performed using the formalism developped in this paper 
(one only needs to include higher multipoles in the sums 
that appear in the expressions of the covariance matrices). 
However, for practical purposes they are probably irrelevant 
for clusters of galaxies. Indeed, higher-order multipoles are 
increasingly difficult to measure and we find that the third 
mutipole, 21 = 4, may only be detected (at a 1 — a level) 
by Euclid or a similar survey. 

We obtain a reasonable agreement with numerical sim- 
ulations for the mean correlations and covariance matrices 
on large scales (r > 10/i _1 Mpc). Redshift-space distortions 
amplify the monopole correlation by about 10 — 20%, de- 
pending on the halo mass. The covariance matrix is also 
amplified by 10 — 40% and the signal-to-noise ratio remains 
of the same order as for the real-space correlation. We find 
that non-Gaussian terms can contribute to the covariance 
matrices up to 20 — 60%, especially for non-diagonal en- 
tries when one of the two bins corresponds to small scales 
(~ 10/i _1 Mpc). They make the correlation matrices less 
diagonal. We also find that the correlation matrices of the 
quadrupole (21 — 2) are significantly more diagonal than 
the correlation matrices of the monopole (2£ = 0). On 
the other hand, the cross-correlations of the monopole and 
quadrupole are rather small, on the order of 0.05. This 
means that the full covariance matrix is approximately 
block-diagonal in the space {Ri,2£i}, that is, we may de- 
couple the analysis of the monopole and quadrupole. 

Our forecasts for on-going and future surveys (only tak- 
ing into account the sample- variance and shot-noise contri- 
butions) show that the impact of redshift distortions on the 
spherical-averaged correlation (i.e. the monopole) will only 
be significant for DES, Erosita, and Euclid. For the other 
surveys investigated here, XXL, SPT, and Planck, the dif- 
ference (Ca o) — (£h«)> between the redshift-space and real- 
space correlations, falls below "1 — <r" of the statistical error 
bar. In practice, even for these surveys it would be better 
to compare the data with theoretical predictions for the 



redshift-space correlations rather than the real-space corre- 
lations, since simple models for the latter arc not more dif- 
ficult to implement and one does not need a very accurate 
modeling of redshift distortions. For covariance matrices, 
it is sufficient to use the real-space ones, which are eas- 
ier to compute. For DES, Erosita, and Euclid, one should 
take into account redshift-space distortions for both the 
mean correlations and their covariance matrices. However, 
because of their modest impact on the covariance matri- 
ces (see Figs. [3] and 0]) and the limited accuracy of such 
computations, it is probably sufficient to use the real-space 
covariance matrices in a first step. 

The quadrupole, 2£ — 2, is more difficult to measure 
from the distribution of clusters of galaxies, and we find 
that only DES, Erosita, and Euclid, among the surveys 
considered here, should provide a useful measure. In this 
respect, clusters are less efficient than galaxies to mea- 
sure redshift-space distortions. Even though clusters have 
a higher bias and larger correlations, they are much less 
numerous than galaxies. This leads to a smaller signal-to- 
noise ratio for the quadrupole. Nevertheless, such a measure 
would provide a nice supplement and would help tighten 
the constraints on cosmology. An extension of this paper 
would be to study on a quantitative level the expected im- 
provement brought by this measure for these surveys. This 
requires a statistical analysis, using Fisher matrices or like- 
lihood functions, that we leave for future works. 

The last multipole generated at linear order by the 
Kaiser effect, 2£ = 4, is very difficult to measure from clus- 
ter surveys and we find that only Euclid should be able to 
provide a meaningful measure. Nonlinear corrections should 
generate higher-order multipoles but they are likely to be 
even more difficult to measure. Therefore, for practical pur- 
poses redshift-space distortions should only impact cluster 
correlations at the level of the monopole and quadrupole 
components. 

We note that our results deliberately disregard obser- 
vational constraints related to cluster redshifts estimation. 
In practice, measuring cluster redshifts is a costly task that 
consists in obtaining spectra for a few (> 10) galaxies per- 
taining to each cluster, leading to errors on the order of 
a(z) < 0.01(1 + 2). This precision depends on the num- 
ber of measured spectra, which in turn depends on cluster 
richness (and consequently, on its mass), cluster apparent 
size, morphological complexity, dynamical state and dis- 
tance to the observer, as well as on fixed instrumental con- 
straints. In some cases, only the brightest central galaxy 
(BCG) can have its redshift measured. This may result in 
a biased cluster redshift estimate, for the BCG may have 
a small proper motion within the cluster potential well 
(lAdami fc Ulmerl 120001: lOegerle fc Hull l200lt ICoziol et al.l 
2009). In many other cases, only photometric redshifts will 
be obtained, resulting in larger reds hift uncertainties - o n 
the order of a(z) ~ 0.05(1 + z) (e.g. iPillepich et alll2012ft . 
However, recent studies have underlined the potential of 
photometric surveys for studying redshift space distortions 
at the galaxy scale, using th e information in t he (2D) angu- 
lar correlation function (e.g. lRoss et al.|[20lTh ■ We leave the 
application of this method to clusters to future work, and 
note that a quantitative analysis (e.g. a Fisher analysis) 
would largely benefit from a proper modelling of follow-up 
incompleteness, uncertainties and biases on cluster redshifts 
given the halo masses and positions. 
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where Si(z) = f dt sin(i)/t is the sine integral. Then, writing the 
two-point correlation function in Eg. 1140 \ in terms of the power spec- 
trum J9}, the angular integration gives 



Appendix A: Computation of geometrical averages 

We give in this Appendix the explicit expressions of the various ge- 
ometrical averages of two-point correlations that are needed to com- 

(s) 

pute the covariance matrix C a - ! , within our framework. We follow the 
method and the notations of Valageas ct al. ( 201l|) as we extend their 
results to redshift-space (especially their Appendix F). 

The second term in Eq. d 3411 involves the average (140 6 of the corre- 
lation between two concentric spherical shells. Defining the 3D Fourier 
transform of the 2£-multipole top-hat as 
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(A.10) 



Next, the last term in Eq J34H involves the non-connected four- 
point average 
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(A.l) 



where again fi T = (e z ■ r)/r, the 3D Fourier-space window of the i-shell 
reads as 
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In Eq JA.ll l the subscript 3 recalls that we consider a 3D window. 

To simplify notations we did not recall this fact in the window 
associated with each 3D radial bin. In studies of the real-space correla- 
tion llValageas et al.ll2011f) . only the usual top-hat window W% (kR) 
appears, which does not depend on the angle of k. Here, because we 
study the anisotropic redshift-space correlation and its various multi- 
poles (55} higher-order multipoles W/g 2 ^(kP) appear, which depend 



(A.3) 



on the direction of k. Eq JA.lll also reads as 
W^ e) (kR) = W^ 2l) (kR) P 2f (Mk), 

where the angular dependence has been factored out, with 

wi 2l \kR) = (-If (M+vJ* j2l{kT) - 
This yields for the window of the i-shell, 
W^ 2e) (k) = W^ 2e) (k)P 2e (^ k ), 



with 
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(A.5) 
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In practice, the redshift-space correlation of clusters of galaxies should 
be dominated by the Kaiser effect and the lower multipoles 21 = 0, 2, 
and 4 as in Eos.l llllH14ll , Therefore, we only give explicit expressions 
for low-order multipoles, which read as 



W?\kR) = 3 Sin(fci?) -^ COs(fcR) , 



W!i 2 \kR) = —^-[kRcos(kR)-4:sin(kR) + 3Si(kR)], 



(A.7) 



(kR) 3 



(A.8) 



MM 



T>i J (Afi) 2 



-(4^ + 1) 



x (Aij +i)p 2l% (av )PM j (h> ) e$e$ • 



(A.H) 



As in|Valagcas ct al. (2011), expressing the two-point correlation func- 
tions in terms of the power spectrum, using the flat-sk y (small an- 
gle) approximation, as well as L imber's approximation (Limberl ll953l : 
Kaiser 1992; Munshi ct al. 2008), we obtain 



e$4'i = % / dk 1 dk 2 p(fc 1 )(i+^) 2 P(fc 2 )(i+/3 M |) 2 



xMfci|| + k n )W^ ) (k 1 )Wp ) (k 2 ) 

x w 2 (\k 1± + k 2± \ve s ) 2 , 



(A.12) 



where fii = (e z • k;)/fc;, 8 a is the angular radius of the survey window 
(in the flat-sky approximation), 

(AO) = tt6I 2 , (A.13) 
and W 2 is the 2D Fourier-space circular window, 



w 2 (k ± m a ) 



do_ iki . De = 2j 1 (k ± ve s ) 



t6» 2 



fei ve a 



(A.14) 



Introducing an auxiliary wavenumber k^ with a Dirac factor 
J'dk^(5£)(kx — kix — k 2 x) and using the exponential representa- 
tion of Dirac functions, we can partially factorize the integrals and 
perform a few angular integrations. This yields 

/ l ( t (r)X ( 2 ' (r) AW (—) . (A. 15) 



MM 

i\3 



2nrV 2 9, 
where we introduced 



/ r \ = rVBs 
\V6sJ ~ 2n 



dk ± e ir -L' k ^ W 2 (kj_V9 s ) 2 



and 



4«,( r )= / dke ik - r P(fc)(l + ^) 2 W-f^)( k ). 



(A.16) 



(A. 17) 



Using the multipole decomposition mil l and the expansion of plane 
waves on spherical harmonics, we can write X?^„ (r) as 

4 



W^\kR) 



27 



- [(105-2(A;P) 2 )fePcos(fci?) 



2(kR) 5 

(105-22(fcP) 2 )sin(fcP) + 15(fcP) 2 Si(fc-R)] , 



(A.9) 



where the sum over I runs up to 4 because the expansion 1111 only 
includes the terms < I < 2 and we only study estimators £)j 2 g 
with < ti < 2. If we consider a model where the redshift-space 
correlation includes multipoles up to 2£ max and we study estimators 
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up to order 2^ m ax the sum (IA. 181) runs up to £ < I 
we introduced 



max i <-max 



Here 



and 



(4m) 



dxx^ 2e ">(x) 



2 



(M) = ^ 

x / duu 2 H/ 3 (24) («) j 2e (au)j 2e ,(bu). 



(A.19) 



(A.20) 



We now turn to the contribution II39II that arises from the four- 
point correlation function, proceeding in a similar fashion. The new 
geometrical average involved by the first two terms within the brackets 
in Eq J39l l reads as 



dr r 



A (t,i')(_^_\ 

VDOJ 



1=0 t'=0 



(A.26) 



which does not depend on the bin i. The third term involves the 
cylindrical average 



The fourth term, which does not factor, gives rise to 



(A.27) 



Then, substituting the expansion HA.18I I into Eq.l lA.15l l we obtain 



drr v i^^wi'^^w 

l,l'=Q 



' (do) ' 



with 



(A.21) 



(A.22) 



where we used that A( s > defined in Eq. IIA.16ll also writes as 

A< s )(y) = A^(y^ 2 )=y I duu W 2 (uf J (yuy/ 1- li 2 ). 
Jo 



(A.23) 



These expressions closely follow those obtained in IValageas et alJ 
1120111 1. To recover their real-space results, it is sufficient to set /3 = 
in Eq.Q, which makes the power spectrum isotropic so that only 
monopole contributions (i.e., £ = 0) remain. Thus, in Eo JA.211) there 
only remains the term £ = £' = 0, and the function A^°'°^(y) can 
actually b e computed explicitly (see their Eqs.(FlO) and (Fll)). 

As in[Valagcas ct al. (2011), we prefer to use the configuration- 
space expression JA. 211 ) rather than the Fourier-space expression 
(IA.12b for numerical computations. Indeed, the integrals have been 
partially factored out in Eo JA.211 1 (so that it has the same level of 
complexity as a two-dimensional integral) and the oscillatory kernels 
jy( 2 4) Q £ Eq. jA.12t have been replaced by the slowly- varying kernels 



W 



of Eo. lfA~T9t . see Eqs. (lB37t - (lR20l) . 



We now turn to the contribution (13S I t that arises from the three- 
point correlation function. The first term within the brackets in 



Eg . (I38t) involves the product of two average s and £\, , which 

are given in Eq.j27j. Following the analysis of Valageas et al. (2011), 
the second term can be written in configuration space as 



^(44 + l)P 2 , i (^ w (r)i;;4(r) 
Using the expansions (1111) and (I A. 18b this yields 



dr 



(A.24) 



dr ; 



t-M c. .A s ) — iq 



x ZV 



E E 



1=0 l'l"=0 



ffW^W^'^), (A.28) 
where we introduced, as in Eo. llA.22t . 



A W) {y) = I J dt,P 2l ^)P 2el (t,)P 2t „^)A^(y,^). 



(A. 29) 



The last two terms within the brackets in Eq. J 39 1) involve the geomet- 
rical average 

X t\%,{v)&\t')&\t + t')A^( < ^\ . (A.30) 

To integrate over angles we again introduce a Dirac factor 
J dr" 5o (r" — r — r') that we write under its exponential representa- 
tion, which we expand over spherical harmonics. Using the expansions 
Ull) and <A.18t we obtain 



7 ,41,11,7 



dr 



R i,+ dr ' 3r '2 /"'' r'' Ar „ r „ 



(™s) 2 A_ R i,+ ~ R i,- J\ r -r'\ rr ' 
E E ^%f{r)&-' 2l 'Hr')&-' U "\r") £ 



1=0 e',t"=o 
n'+e" 



E (-l) n+n ' + ^"(4£ l + l)(4n+l)(4n' + l) 



n— I n 

' 2n 



' 24 2£' V / 2n In' 21" V 4 (f.n) f r \ 
; V \-dbJ 



X C n ' n '^"(r,r',r"), 



where we introduced the symmetric functions 



(A.31) 



As) As) 



«»+i) 



dr 3r 2 



i? 3 



R 3_ 

8 ' f=0 «'=0 



Hf^ ( 



2£ 2£' 24 




(A.25) 



The third term within the brackets in Eg. (138D is obtained from 
Eo JA.25ll by exchanging the labels "i" and "j". 



C«i,«2,«3( ri)r . 2!r . 3 ) = - ri r 2 r 3 / dk k 2 j 2ei (kr!) 

Jo 

x j 2 i 2 {kr 2 )j 2l Jkr s ). (A.32) 



For numerical computations it is more efficient to write Eg. (IA.31I I in 
a partially factorized form by exchanging the order of the integrations 
and moving to the left the integration over k included in the function 
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C*i'*2i*3. Using the decomposition (1 1 1 ft - <l 1411 and integrating over r" 
this gives 



A") A s ) As) 
j f ;i i' ;j 



20.. 



4 2 li+f n'+t" 

t^E E E E 

1=0 ti,t"=On'=\li-l'\n=\n'-t"\ 
x(-ir +n ' +< "(« l + l)(4n+l)(4n' + l)( 2 ^' 2 J ^ J 

, o poo 

N ' dr r ■ lu \T^ 2i )f \ 



( 2n 2n< 2t" \ 2 m » 
\ o J 



(A.33) 



where F( 2t ) is the prefactor in Eg s, 1121- 11141 

2 4/3 4/3 2 



p(0) = i + M + <L F W 
3 5 ' 



F (4) = 2^±_ 

7 ' 35 



Again, all these geometrical averages of redshi f t-spac e correlations 
reduce to the real-space results of Valageas ct al. (2011) when we set 
/3 = and only keep the monopole terms. 



Appendix B: Some useful integrals 

We describe here the computation of some of the functions introduced 
in Appendix lAl 



B.l. Multipoles ofA^iy) 

Using the expression HA.23I I and the expansion ( Gradshteyn & Ryzhikl 
1991) 



we can write 



^(2^ + 1/2) 



(2n - 1)!! 
2"n! 



x 4+l/j(z)ftt.M. 



(B.l) 



B.2. Symmetric functions C l1,12,13 

The functions C^ 1 '^ 2 ^ 3 defined in Eo. llA.32ll vanish if the three lengths 
{ri, r 2 , r$\ do not obey triangular inequalities. It is also useful to note 
that these functions obey the scale invariancc 

A>0: C^^ 2 ^ 3 (Ari,Ar 2 , Ar 3 ) = C fl ^ 2 ' f3 (ri,r 2 ,r 3 ). (B.7) 

For |ri — r 2 | < rz < r% + r- 2 , as in Ea. llA.31l l. we obtain explicit 
expressio ns using the recursion relatio ns of Bessel functions and the 
property l IGradshtevn & Ryzhik 19941) 

(bc) a ~ 1 

dkk 1 - a J a {ak)Jp{bk)J l3 {ck) = - J_ - (sin«) a ~ 1/2 



(B.8) 
(B.9) 



x J2i 2+ i/2(.kr 2 ) Ju 3 +l/2{krz). (B.10) 
Then, using the recursion 



2na° 
X P p-l/2 ( cosl, )> 

for 

\a - b\ < c < a + b, 2bc cos(v) = b 2 + c 2 - a 2 . 
Indeed, Eq.l lA.32ll also reads as 



(A.34) o«i,< a A( Pl 

, r2,rz) = \/2-Kr\r2r-i I dk Vk J 2 e 1+1 /2(kri) 



4f — 1 

^+l/2( 2 ) = — " — J2l-l/2( z ) - J2t-3/2( z ), 



(B.ll) 



we can lower the indices {^1,^2,^3} and we can express all required 
quantities in terms of integrals of the form (IB . 8l t with two identical 
indices and the right power-law prefactor. For instance, the first few 
functions read as 



c u,u,o = x 



C 



o,r,i _ 1 



C 



1,1,1 _ 3 r 2 r 3 



2 ri 



3ti 2 - 1 



were wc defined 



Mi 



r 2 + r 3 ~ r l 
2T2T3 



(B.12) 
(B.13) 

(B.14) 
(B.15) 



A^Hv, ^) = ^(4n+l) (2 ^^ )!! B n (y)P 2n (»), (B.2) 

n = 

with 



Bn(y) = V I duuW 2 {u) j2n(yu). 

'o 



This yields the multipoles 

= V" (4n+l) 

n=|i 1 -^ 2 | 

x B n (y), 



(2n-l)!! / 2h 2i 2 2n V 
V J 



2"n\ \ 



and 



^1+^2+^3 



(B.3) 



(B.4) 



A (h,l2,e 3 ) {y)= (4n+l) (2?1 1)!! y^^ 2 ' f3 -"B n fa), (B.5) 

Z — / 2 n ri! 



with 



) 1 ' ~ J '■" = I / d M P 2 ^ (n.)P« 2 (^p^ ( M )P^ 4 (B.6) 

Because we only consider the low-order multipoles < £i < 2 of the 
redshift-space correlation function, we only need the first few mul- 
tipoles A(*i^2) and A^i^ 2 '^ 3 ). The associated numbers Y^M^zM 
and functions B n (y) can be computed in advance. 



B.3. Functions W/ (a, 6) 

To compute the geometrical means I. (r) of the redshift-space cor- 
relation with respect to bin i, defined in Eg. (IA. 17t . we need the func- 
tions W^' defined in Eq JA.2011 . Using Ea JA.4l l. they also write as 



W. 



f')(a, b ) = (-1)^(44 + 1) * 



dr — C^'** (a, b,r), 



(B.16) 



where the functions C^ 1 '* 2 '^ 3 were defined in Eq. IIA.32ll . From the 
triangular constraint on these functions C 1 * 1 '^ 2 '^ 3 we obtain 

[o-6[ >1: Wff ] =0, (B.17) 
while the scale invariance JB. 711 leads to 

a + b<l : Wf^'\a,b) = {a + b) 2 wf^ f_2_,_5_V 
1 1 \a + 6 a + b J 

(B.18) 

Then, from the explicit expressions obtained in Appendix IB. 21 which 
are rational functions of {ri , r 2 , r$}, we can perform the integral 
IB. 161 1 and derive explicit expressions. In particular, the first func- 
tion reads as 



|o-6| < 1, a+b < 1 : W { 



(0,0) 



3ab, 



(B.19) 



\a-b\ <1, a + b>l: W^°' 0) = -[1 - (a - b) 2 ]. (B.20) 
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k'°], 2(,=2(,=2, M>2xl0'V ] A<t, 



ft'-j, 2ti=0, 2t/=2, M>2xl0 L V'yVt o 
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the redshift-space auto- 
for the 



Fig. C.l. Contour plots for 

correlation matrices TZ\ S j, defined as in Eq. (j4"Tj) . tor 
u 2l — 2" redshift-space correlation studied in Sect. 14.31 
We consider the same cases (distance bins, redshift inter- 
val, and mass thresholds) as in Fig. [SI where we plotted the 
auto-correlation matrices of the "monopole" term of the 
redshift-space correlation. 



Appendix C: Correlation matrices of the "2£ = 2" 
redshift-space correlation 

C.l. Auto-correlation of the "2£ = 2" multipole 



We show in Fig. IC.ll the auto-correlation matrices of the "2£ = 2" 
redshift-space correlation. For simplicity we only plot our results for 
the full matrices that include the non-Gaussian contributions associ- 
ated with the 3-pt and 4-pt correlation functions. We can see that 
they are significantly more diagonal than the auto-correlation matri- 
ces of the "monopole term" that were displayed in the right panels of 
Fig. El 



C.2. Cross-correlation of the "2£ = 0" and "21 = 2' 
multipoles 



We show in Fig. IC.2I the cross-correlation matrices of the "2£i = 0" 
and "2£j = 2" redshift-space correlations. Thus, in the plane {Ri, Rj} 
we plot the contour plots of the matrix 



7v 



K it- 



c 



o 



,<7 



« 



(C.l) 



Because li ^ £j this matrix is no longer symmetric in the plane 
{Ri,Rj} and it can be negative for Ri = Rj. We can see that its 



plots for the redshift-space cross- 
[s) defined as in Ea. lfCTTj) . for the 



Fig. C.2. Contour 
correlation matrices TZ 
u 2£i — 0" and u 2£j — 2" redshift-space correlations. We 
consider the same cases (distance bins, redshift interval, 
and mass thresholds) as in Figs. [51 and lC.ll 



amplitude is rather small, typically on the order of 0.05, even along 
the diagonal. Indeed, because l{ ^ £j the pure shot-noise contribu- 
tion to Eq. l] .' > 4 t i vanishes (the factor Sg. ^ . (4£; + 1) within the brackets) . 
This is due to the orthogonality of the Legendre polynomials, which 
are the basis of the multipole decomposition of the redshift-space cor- 
relation. The small amplitude of {Ri,Rj} found in Fie. IC.21 means 

that the full covariance matrix ' is approximately block-diagonal 

in the space {Ri,2£i}. Thus, we may neglect ' for ti 7^ £j 

(s) (s) 
decouple the analysis of £ R Q and £ R ' „ . 



and 



References 

Abell, G. O. 1958, ApJS, 3, 211 

Adami, C. & Ulmer, M. P. 2000, A&A, 361, 13 

Benjamin, J., Heymans, C, Semboloni, E., et al. 2007, MNRAS, 381, 
702 

Bernardeau, F., Colombi, S., Gaztanaga, E., & Scoccimarro, R. 2002a, 

Phys. Rep., 367, 1 
Bernardeau, F., Mellier, Y., & van Waerbeke, L. 2002b, A&A, 389, 

L28 

Bernstein, G. M. 1994, Astrophys. J., 424, 569 

Blake, C, Davis, T., Poole, G. B., et al. 2011a, MNRAS, 415, 2892 
Blake, C, Glazebrook, K., Davis, T. M., et al. 2011b, MNRAS, 418, 
1725 

Bohringer, H., Schuecker, P., Guzzo, L., et al. 2001, A&A, 369, 826 
Borgani, S., Rosati, P., Tozzi, P., et al. 2001, ApJ, 561, 13 
Burenin, R. A., Vikhlinin, A., Hornstrup, A., et al. 2007, ApJS, 172, 
561 



19 



P. Valageas and N. Clerc: Redshift-space correlation functions in large galaxy cluster surveys 



Cabre, A. & Gaztanaga, E. 2009, MNRAS, 396, 1119 
Clcrc, N., Sadibekova, T., Pierre, M., et al. 2012, MNRAS, in press 
Cole, S., Percival, W. J., Peacock, J. A., ct al. 2005, MNRAS, 362, 
505 

Colombi, S., Dodelson, S., & Widrow, L. M. 1996, ApJ, 458, 1 
Coziol, R., Andcrnach, H., Caretta, C. A., Alamo-Martinez, K. A., & 

Tago, E. 2009, AJ, 137, 4795 
Dcsjacques, V., Crocce, M., Scoccimarro, R., & Sheth, R. K. 2010, 

Phys. Rev. D, 82, 103529 
Eisenstein, D. J. & Hut, P. 1998, ApJ, 498, 137 

Eisenstein, D. J., Weinberg, D. H., Agol, E., et al. 2011, AJ, 142, 72 
Eisenstein, D. J., Zehavi, I., Hogg, D. W., et al. 2005, ApJ, 633, 560 
Fry, J. N. 1984, ApJ, 277, L5 

Gradshtcyn, I. S. & Ryzhik, I. M. 1994, Table of integrals, series and 
products, cd. Gradshtcyn, I. S. feamp; Ryzhik, I. M. (Academic 
Press, New York) 

Groth, E. J. & Peebles, P. J. E. 1977, Astrophys. J., 217, 385 

Guzzo, L., Pierleoni, M., Meneux, B., et al. 2008, Nature, 451, 541 

Hamilton, A. J. S. 1992, ApJ, 385, L5 

Jackson, J. C. 1972, MNRAS, 156, IP 

Kaiser, N. 1987, Mon. Not. R. Astron. Soc, 227, 1 

Kaiser, N. 1992, ApJ, 388, 272 

Katgcrt, P., Mazure, A., Perea, J., et al. 1996, A&A, 310, 8 
Kilbinger, M., Benabed, K., Guy, J., et al. 2009, A&A, 497, 677 
Koester, B. P., McKay, T. A., Annis, J., et al. 2007a, ApJ, 660, 239 
Koester, B. P., McKay, T. A., Annis, J., et al. 2007b, ApJ, 660, 221 
Komatsu, E., Smith, K. M., Dunkley, J., et al. 2011, ApJS, 192, 18 
Landy, S. D. & Szalay, A. S. 1993, ApJ, 412, 64 
Limber, D. N. 1953, ApJ, 117, 134 
Manera, M. & Gaztanaga, E. 2010, arXiv:0912.0446 
Mantz, A., Allen, S. W., Rapetti, D., & Ebeling, H. 2010, MNRAS, 
406, 1759 

Massey, R., Rhodes, J., Leauthaud, A., et al. 2007, ApJS, 172, 239 
Mehrtens, N., Romer, A. K., Lloyd-Davies, E. J., et al. 2011, ArXiv 
e-prints 

Meiksin, A. & White, M. 1999, Mon. Not. R. Astron. Soc, 308, 1179 
Melin, J.-B., Bartlett, J. G., & Delabrouille, J. 2006, Astron. 

Astrophys., 459, 341 
Moles, M., Bem'tez, N., Agucrri, J. A. L., ct al. 2008, AJ, 136, 1325 
Munshi, D., Valageas, P., van Waerbeke, L., & Heavens, A. 2008, 

Phys. Rep., 462, 67 
Nishimichi, T. & Taruya, A. 2011, Phys. Rev. D, 84, 043526 
Oegerle, W. R. & Hill, J. M. 2001, AJ, 122, 2858 
Pacaud, F., Pierre, M., Adami, C, ct al. 2007, MNRAS, 382, 1289 
Pacaud, F., Pierre, M., Refregier, A., et al. 2006, MNRAS, 372, 578 
Peacock, J. A., Cole, S., Norberg, P., ct al. 2001, Nature, 410, 169 
Peebles, P. J. E. 1980, The large-scale structure of the universe, ed. 

Peebles, P. J. E. (Princeton University Press, Princeton, N.J., USA) 
Peebles, P. J. E. & Hauser, M. G. 1974, Astrophys. J. Supp., 28, 19 
Percival, W. J. & White, M. 2009, Mon. Not. R. Astron, 393, 297 
Pierre, M., Pacaud, F., Juin, J. B., et al. 2011, MNRAS, 414, 1732 
Pillepich, A., Porciani, C, & Reiprich, T. H. 2012, MNRAS, 422, 44 
Planck Collaboration, Ade, P. A. R., Aghanim, N., et al. 2011, A&A, 

536, A8 

Predehl, P., Boehringer, H., & et al., H. B. 2009, Proceedings of the 

conference "X-ray Astronomy 2009", Bologna, September 2009 
Refregier, A., Amara, A., Kitching, T. D., et al. 2010 
Reichardt, C. L., Stalder, B., Bleem, L. E., et al. 2012, ArXiv e-prints 
Reid, B. A., Samushia, L., White, M., et al. 2012, ArXiv c-prints 
Reid, B. A. & White, M. 2011, MNRAS, 417, 1913 
Ross, A. J., Percival, W. J., Crocce, M., Cabrc, A., & Gaztanaga, E. 

2011, MNRAS, 415, 2193 
Ross, N. P., da Angela, J., Shanks, T., et al. 2007, MNRAS, 381, 573 
Rozo, E., Wechsler, R. H., Rykoff, E. S., et al. 2010, ApJ, 708, 645 
Sarazin, C. L. 1988, X-ray emission from clusters of galaxies 
Schrabback, T., Hartlap, J., Joachimi, B., et al. 2010, A&A, 516, A63 
Sehgal, N., Trac, H., Acquaviva, V., et al. 2011, ApJ, 732, 44 
Slosar, A., Font-Ribera, A., Pieri, M. M., et al. 2011, J. Cosmology 

Astropart. Phys., 9, 1 
Smith, R. E., Peacock, J. A., Jenkins, A., et al. 2003, MNRAS, 341, 

1311 

Spergel, D. N., Bean, R., Dore, O., ct al. 2007, ApJS, 170, 377 
Sunyacv, R. A. & Zeldovich, Y. B. 1972, Comments on Astrophysics 

and Space Physics, 4, 173 
Szapudi, I. 2001, in Annals of the New York Academy of Sciences, 

Vol. 927, The Onset of Nonlinearity in Cosmology, ed. J. N. Fry, 

J. R. Buchler, & H. Kandrup, 94 
Szapudi, 1. & Colombi, S. 1996, ApJ, 470, 131 



Taruya, A., Nishimichi, T., & Saito, S. 2010, Phys. Rev. D, 82, 063522 
Teyssier, R. 2002, Astron. Astrophys., 385, 337 
Teyssier, R., Pires, S., Prunet, S., et al. 2009, A&A, 497, 335 
Tinker, J., Kravtsov, A. V., Klypin, A., et al. 2008, ApJ, 688, 709 
Tinker, J. L., Robertson, B. E., Kravtsov, A. V., et al. 2010, ApJ, 
724, 878 

Valageas, P., Clerc, N., Pacaud, F., & Pierre, M. 2011, A&A, 536, 
A95 

Vanderlinde, K., Crawford, T. M., de Haan, T., et al. 2010, ApJ, 722, 
1180 



20 



